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CN ! Abstract 

The signal strength of the recently discovered Higgs boson-like particle in the 
■ diphoton channel seemingly constrains physics beyond the standard model to 

a severe degree. However, the reported signal strength is prone to possible un- 
derestimation of uncertainties. We propose a discriminant that is relatively free 
of many of the theoretical uncertainties, and use this to gauge the impact on 
the phenomenological MSSM. A Bayesian global fit to all the pre-LHC data 
results in posterior distributions for the masses that are neither very restrictive, 
nor sufficiently prior-independent (except for the Higgs and stop masses). The 
p I imposition of the Higgs data, on the other hand, yields interesting and nearly 

' . prior-independent constraints. In particular, the existence of some light super- 

partners is favoured. 

1. Introduction 

The recent discovery of a 126.5 GeV or 125.0 GeV Higgs boson-like state, respectively by the ATLAS 
^ and CMS collaborations lUEl, at the LHC marks the beginning of an exciting scientific era for under- 

^ . standing the mechanism responsible for electroweak symmetry breaking and for endeavour in seaixh for 

\ new, beyond the standard model (BSM), physics. For a standard model-like Higgs boson, the search 

O ■ strategies at the LHC involve, in order of decreasing sensitivity, production processes driven by gluon 

fusion gg — )• H, vector boson fusion qq — )• Hqq, associated production with a vector boson qq — )• VH 
and the associated production with a top-quark pair gg — )• tlH. Further to production, the main seaixh 
channels are defined by decay modes into 77, ZZ* , WW*, bb and r+r^. Subsequent decays of the W- 
and Z-bosons with at least one leptonic decay are also considered. For the experiments with the Tevatron 
collider, the main production modes, for the most sensitive search topologies, are the associated produc- 
ed I tion of a Higgs with a vector boson qq — )• VH, with the subsequent decay of the Higgs into a bb pair, and 
the fl'g-fusion production with the subsequent decay into two W bosons which then decay leptonically. 
The LHC discovery has supporting evidence from the Tevatron experiments (see e.g. |f3l and references 
therein) where the excess seen is in accord with a 115 to 135 GeV Higgs boson. 

Usually, the Higgs boson signal is quantified by the production cross section times the relevant 
decay-channel branching fraction to give the signal strength, 

ajgg ^ h) Brjh ^ X) 
a{gg h)sM Br{h X)sm ' 

Here X = 77, ZZ* , WW*, bb, or t~^t~. Of the Higgs boson decay channels reported by the experi- 
mental collaborations, fi^^ apparently makes for a very strong BSM-constraining observable. However, 
as argued in Ref. H, it is relevant to question whether the apparent excess is actually due to new physics 
effect or due to under-estimations in theoretical uncertainties. The question would also propagate di- 
rectly to any analysis that employs the reported signal strength. In order to avoid this and for the purpose 
of allowing robust conclusions to be derived from the discovery data, we ai^e going to propose a discrim- 
inant between new physics and the standard model (SM), which should be largely free of the mentioned 
source of uncertainties. 
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We consider, here, the R-parity conserving minimal supersymmetric standard model (MSSM) as 
a template for BSM. Rather than appeal to a particular scenario for supersymmetry breaking and/or 
make an Ansatz for the parameters at some high scale, we, instead, use the phenomenological MSSM 
(pMSSM) parametrisation ||5l,l6l to characterise the supersymmetric sector. Although this is almost the 
most general scenario given the filed content, it is, nonetheless, subject to very many constraints. For 
example, Refs. |17J[8l had used all pre-LHC data to aiTive at Bayesian posterior distributions for the 
superpartner masses. We discuss here how the mass of the Higgs resonance and the information on the 
signal strengths (in particular, the robust discriminant) can be used to further refine the posterior to yield 
some very interesting clues about the pMSSM spectrum. Other discussions concerning LHC constraints 
on the pMSSM were effected in Refs. |[9UT3l. 

Several studies have been made of the MSSM Higgs sector in light of the observed excess (Refs. (13\ - 
\VJ\ comprise a non exhaustive list), in the diphoton channel signal strength or the effect of light sparticles 
in altering it. Some of these attempts are based on random scans in the pMSSM parameter space while 
others restrict themselves to sub-spaces or hyper-planes of the MSSM parameters space. In this article 
we make use of a complete set of pMSSM points that came out of statistically convergent Bayesian fits. 
The aim is to derive the implications of the recent LHC discovery on the pMSSM posterior distributions 
by gauging its effect on the sparticle masses. In the the rest of this section we introduce the BSM discrim- 
inant for the analysis. In Section|2]] we give a brief recapitulation of the pMSSM's posterior distributions 
for sparticles masses, Higgs boson signal strengths, and the proposed discriminant predictions from a 
pre-LHC Bayesian fit to data. Section [3] describes the implementations and implications of the Higgs 
boson data on the pMSSM predictions. 

1.1 A useful discriminant 

The ratios fix, see eq.©, represent good measures for possible deviations from the SM and have, indeed, 
been exploited to discriminate against BSM models. However, such an exercise ought to be approached 
with care, for it is yet early days of Higgs physics and the statistics are still low. A naive adoption of the 
current bounds would indicate a significant tension with the SM itself ifTSl . Of particular concern are the 
theoretical eiTors pertaining to the choice of the paiton distributions and the renormalization (Qr) and 
factorization (Qf) scales. These uncertainties can be large and have been addressed in Ref. |[T9l . The 
LHC Higgs cross section working group recommends EOl that in convoluting the errors, these should be 
added linearly. However, both ATLAS and CMS have added them in quadrature. As Ref. H explicitly 
demonstrates, adopting the more conservative approach of linear addition would substantially increase 
the eiTor bars to the extent of a la agreement with the SM. A further uncertainty pertains to the use 
of an effective field theory (in the Umit of an infinitely heavy top) used for many of the higher-order 
calculations. 

This criticism brings into question the naive use of fix- However, note that the uncertainty applies 
primarily to the production process under consideration viz. pp ^ h (driven, at the lowest order, through 
gluon-fusion). In fact, as far as non-hadronic final states (diphoton, and four-lepton, charged or neutral) 
are concerned, this is an exact statement as quantum chromodynamics (QCD) coiTcctions can occur 
exclusively in the initial state alone. For — )• /i — )• bb, though, one could have gluon exchanges 
between the initial and final state particles, thereby engendering non-universal dependence on the parton 
density. Nevertheless, as is well-known, such diagrams result in only a subdominant component of the 
higher-order QCD corrections, and, thus, even for the pp ^ h ^ bb case, the dependence is identical 
to a very good approximation. The same ai^gument holds for pp ^ h ^ VV* with the gauge boson(s) 
going into jets. 

In other words, if we consider only those final states that are primarily driven by Higgs-production 
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Search channel 


ATLAS /i at m/, = 126.5 GeV 


CMS /i at m/, = 125.0 GeV 


i/ — > 77 
H ZZ* 


1.8 ±0.5, 4.5c7 
1.2 ±0.6, 3.6cT 


1.6 ±0.4, 4.6cr 
0.71° j, 3.2a 



Table 1: The ATLAS m and CMS IS magnitudes for the Higgs boson production signal strengths. 



through gluon-gluon fusion, the ratios of ratio^, namely, 

Rij = — , (2) 

are essentially free of the aforementioned uncertainties due to the choice of Qr^f as well as the choice 
of parton distribution^ Here i and j represent distinct choice from X = 77, ZZ* , WW*, bb, or r+r". 
The same argument would hold for other systematic errors that are factorisable, or are nearly universal. 
A prime example is afforded by the uncertainty in the luminosity measurement. 

Thus, it makes eminent sense to consider all the independent ratios Rij, and, at best, one individual 
/ij, for, together, they embody the same information as the set {/Uj}. There is, of course, a price to pay: 
the relative statistical eiTors in Rij would be larger than those in the individual fii. Furthermore, certain 
statistical eiTors in /^tj are con^elatecH, and, unless the coiTclation matrix is available, a naive combination 
would tend to overestimate the error in Rij. This, certainly, is the case in the current context in the 
absence of sufficient pubhc information about the ATLAS and CMS analyses. 

Nonetheless, for our analyses we shall consider such ratios Rij. Of particular interest here is the 
ratio of the two reported discovery channels, namely four-lepton and diphoton. The ATLAS and CMS 
ratios for the same, as shown in Tabled] are consistent with each other within about 0.8a, and combining 
the two, we have, 

i?^^.^^ :^ ^ = 0.56 ±0.25 , (3) 

which, in itself, is consistent with the SM within 2a. Yet, the very fact that the central value of Rzz ;77 
is substantially below the SM expectation would turn out to have rather interesting consequences. It 
should be realized that this deviation is primarily driven by the fact that ^u^^ is significantly larger than 
unity; yet, unlike the latter, this observable is largely free from the aforementioned theoretical errors. 
Moreover, the enw in eq.© has been somewhat overestimated and would be reduced substantially once 
the eiTor correlation matrix is available. This renders any inference drawn from the use of eq.© to be 
conservative in nature. 



Similarly, we also have 



Rzz;ww^^^ = 0-95lt^oill, (4) 



which is rather consistent with the SM. And while pp ^ H ^ t~^t seems to be much smaller than the 
SM expectations, the current eiTor bars are too lai^ge for the observable to be of any consequence. 

2. A brief recapitulation of the pMSSM and its predictions 

Within the MSSM, the masses and interactions of the Higgs bosons are strongly affected by loop cor- 
rections, in particular from the sector of the third-generation quarks/squarks. Here, we consider the 



'Ref. 1211 too had raised the possibility of using such ratios, without putting it into practice. 

^It is also worthwhile to note that even the ratios of, say, gluon-fusion and vector-boson fusion cross-sections also profit on 
this account, though not to the same high degree as is applicable here. 

''An example would be those pertaining to the electromagnetic calorimeter measurements, relevant to both the diphoton 
channel and the VV* modes containing electrons. 
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predictions for the Higgs bosons in the frame of the phenomenological MSSM (pMSSM) where a non 
biased approach to the MSSM paiametrisation is adopted and with the parameters defined and simultane- 
ously varied at the electroweak scale, the energy scale which the LHC is now probing. For the pMSSM 
fit, the parametrisation is completely decoupled from the details of the physics responsible for the break- 
ing of SUSY. Only real soft SUSY breaking terms were considered, with all off-diagonal elements in the 
sfermion mass terms and trilineai^ couplings set to zero, and the first-and second-generation soft terms 
equalisecill leading to a set of 20 parameters: 

M ■ ^Srdgen lst/2ndgen_ . 9 , „^ o 

}q,u.d,l.e JQ,U,D,L,E 

where Mi, M2 and M3 ai^e the gaugino mass parameters; mj are the sfermion mass parameters. f, ,-^^=e 
represent the tiilinear scalar couplings while the Higgs-sector parameters are specified by , rn'jj^ , 
tan/3 and sign{/j,). The set of observables used for the pre-LHC global fits ITHH are the CDM relic 
density, electroweak physics observables (EWPOs) and B-physics observables: 

O = {mw, sin2 elj^, Tz, 6a^, Rf, A%, £ = A', R^^, A)1 A^'^ (6) 
BR{B X,7), BR{Bs Aq-, Rbr{Bu^tu), RaMb,, 

^CDMh^}- 

In the absence of direct data, the pMSSM parameter space is weakly constrained in the sense 
that only a few of the posterior distributions (in particulars those that are strongly linked to the Higgs 
sector) are approximately prior-independent. The Bayesian fits to data, nonetheless, have an important 
virtue of indicating how effective a particular data is in constraining a given model in an unambiguous 
manner. For example, if the LHC data rules out a prior-independent pMSSM prediction, then a large 
class of MSSM possibilities would have been ruled out at a go instead of trying many little constructed 
sub-MSSMs. This is what we set to explore in our analysis: (i) check if the Higgs boson discovery data 
rules out the approximately prior-independent pMSSM predictions in the Higgs sector and (ii) search for 
any prior-independent impact of the data on the pre-LHC posterior samples. 

Unlike the case for sparticle spectrum, the Higgs sector of the pMSSM is strongly constrained by 
the LEP lower limit and the MSSM's theoretical upper bound on the Higgs mass. This is reflected in 
the results for the pre-LHC Bayesian fits ITKH where the Higgs sector results were found to be approxi- 
mately prior-independent. The light CP-even Higgs boson mass was predicted to be between 117 GeV to 
129 GeV at 95% Bayesian confidence level for both logarithmic and flat priors (i.e., approximately prior 
independent). This is consistent with the recent observation of a new state compatible with the Higgs 
boson hypothesis near 125 GeV. The constraint on the Higgs sector, in turn, puts a severe constraint 
on the mass of the scalar top-quark which plays a crucial role in the radiative contribution to the Higgs 
boson mass. As a result, the statistically convergent Bayesian fits to indirect collider and cold dark matter 
relic density data for the pMSSM ||71[8l reveal that the scalar top-quarks masses have an approximately 
prior-independent posterior distribution centred around 2 TeV. 

The Higgs boson mass in the pre-LHC log prior Bayesian fit has a maximal correlation, q, with 
the top-quark mass; followed by the EWPO Rc, the Higgs sector bilinear parameter fi, Q = ^m^^ m^^ , 
^Zrdgen^ TTi^^, uij^'^^'^"', TTT-j^, and BR(6 — ^ S7) with Correlation strengths of 0.71, 0.65, 0.53, 0.49, 0.45, 
0.42, 0.42 and 0.41 relative to q respectively. For the flat prior, case the maximal correlation is with 
the /i pai^ameter; followed by the top-quai^k mass, Rc, tan (3, BR(fe — )• sj), Q = ^m^-^ m^-^, rr?^'^^^"'. 



Tz, mi^, rnj'^^'''^ and m^-^ with correlation coefficients of 0.86, 0.79, 0.65, 0.57, 0.49, 0.46, 0.44, 0.44, 
0.41 and 0.38 relative to q respectively. Unlike the case for the Higgs boson mass, the conformity of the 
pMSSM prediction with the reported signal strengths is not immediately visible. 



'^This set of assumptions are obviously designed to conform to the very strict constraints on CP violation as well as flavour 
changing neutral current processes, whether in meson mixings or in rare decays. 
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We assume that the peak discovered by the ATLAS and CMS collaborations is explained by the 
pMSSM's light CP-even Higgs boson. We then compute the signal strengths for the production of the 
Higgs boson ail 8 TeV LHC. These are used to check the compatibility and implications of the discovery 
data on the pMSSM's both approximately prior-independent (such as for the scalar top-quarks) and prior- 
dependent (as is the case for other sparticle masses) posterior predictions. 

2.1 The pre-LHC predictions for ^ix and Rzz m 

For each pMSSM model, the Higgs boson production cross section a{g g h) and their branching 
fractions Br{h — X) are computed using the high energy physics package FeynHiggs |[22l . From 
these, the signal strengths /ix eq-ffl are calculated. The pMSSM posteriors for the signal strengths are 
shown in Fig.[T] The signal strengths have a spread encompassing unity, but typically have somewhat 
reduced values. In particular, the signal strengths for gg ^ h ^ bh/r^T^ both, lai^gely, tend to be below 
those expected for the standard model. Naively, this might seem unexpected in view of the fact that the 
Yukawa coupling in the MSSM is often larger than in the SM (at least for large tan /3). However, this has 
to be convoluted with the fact that the gg ^ h amplitude is suppressed on account of the cancellation 
between the top and the stop loops. While this cancellation is less effective for super-heavy stops, such 
solutions are, typically, less favoured once constraints such as ^cDMh"^, BR{B — Xs^) and are 
taken into account. Indeed, this could be gleaned from the pre-LHC pMSSM posteriors as well. The 
reduction is more pronounced for the log prior fit's posterior distribution since the spaiticle masses are 
lighter. Note that — )• /i — )• 66 would, perhaps, never be seen over the large QCD background, while 
both ATLAS and CMS see fewer candidates for — /i — > t^t~ than would be expected within the 
SM. Thus, no conclusion can be drawn, as yet, from these two modes. 

The ZZ mode strength can be quite consistent with the SM value, with relatively fewer points 
leading to a suppression. On the other hand, the observed rate too is marginally smaller than that in 
SM with an error bar sufficiently large to preclude any significant constrain^. The situation with the 
diphoton rate is more dramatic. While, on the face of it, even //^^ can be consistent with unity, note 
that this rate largely lies well below the SM one. Part of this can be attributed to the aforementioned 
suppression in the production cross section. An additional suppression accrues due to the cancellations 
between the VF-loop (dominant contribution within the SM) and the chargino loops (care must be taken 
to include the charged-Higgs loops as well) on one hand and the top- and stop-loops on the other. More 
importantly, the observed rate is significantly above the SM expectation, a situation that is rare within 
the MSSM. Thus, if we use /x^^ as a constraint, very few points would survive in the pMSSM scan, a 
situation that is con^oborated by Table |2j where we have summarised the relative number of the model 
points that survive the various experimental observations, individually and in combination. 

3. Higgs discovery impact on pMSSM sparticle predictions 

The Higgs discovery data may be considered to consist of two forms of observables, namely the Higgs 
boson mass and the signal strengths. Although the pre-LHC data predict m/j G [117, 129] GeV at 95% 
C.L. in correspondence with large stop masses as every additional increment in mh beyond mh = mz 
(with only the MSSM particle content) needs progressively larger stop masses Q This immediately 
implies that requiring m/j = 125 — 126.5 GeV would impose significant further constraints on the pre- 
LHC pMSSM posterio distribution, see Table|2l However, given the about 3 GeV errors from theoretical 
computations |23] one could, for the purpose of the present analysis, consider the Higgs boson mass 

^Although the discovery claims of both ATLAS and CMS are based on a combination of 7 TeV and 8 TeV data, the latter 
has a somewhat weight on account of both higher cross sections as well as a larger integrated luminosity. Furthermore, the 
difference between the two runs accounts for little in our analysis. 

*It should be noted here that the situation with h WW is similar, albeit with even larger error bars. 

'There is an exception for the "Natural SUSY" scenario which allows for light stops and at the same time satisfying 
rrih ~ 125 GeV. This is possible within the pMSSM but is beyond our purpose here. 
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Fig. 1: 2-D marginalised posterior distributions for the log prior pMSSM's Higgs boson signal strengths 
at the LHC with 8 TeV collision energy. The colour spread from light-pink to blue represents relative 
probability spread from zero to unity. The inner and outer contours represent 68% and 95% probability 
regions respectively. The distributions for the flat prior sample is approximately the same but with 
nan^ower peaks around unity. The 1-D marginalised posterior distribution are for the log 

(red/broken curve) and flat (black/solid curve) prior pMSSM fits. 



range of 

rrih = 122 - 128 GeV (7) 

to be plausible and allowed for the pMSSM predictions. 

The experimental values for the Higgs boson signal strengths we use are shown in Tabled] While 
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the excesses do depend on the exact mass bin, it is, nonetheless, useful to combine the two data from the 
two experiments to yield 

1.2<^^^<2.3, 0.4 < < 1.8. (8) 

While this (along with eq.©) does represent a slight simplification, this is the best that can be done until 
the two experiments combine their data. After applying various constraints including the new discovery 
mode signal strengths, the relative amount of the model points that survive various combinations of the 
experimental cuts eq.Q and eq.© are summarised in Table |2l 



No. 


Constraint 


Log prior survive 


Flat prior survive 


1. 


ruh = 122 - 128 GeV 


36.08% 


41.24% 


2. 


nih = 125 - 126.5 GeV 


9.17% 


9.93% 


3. 


Rzz;yy = 0.56 ± 0.25 


0.30% 


0.00% (3 points) 


4. 


2-0" RzZ]-y-y 


99.90% 


100.00% 


5. 


1.61 -a Rzz;-y-r 


13.30% 


0.40% 


6. 


mh{l) & 1.67-0- Rzz-a-y 


4.67% 


0.20% 


7. 


1.2 < /i^^ < 2.0 


0.20% 


0.01% (5 points) 


8. 


0.4 < ^izz < 1-8 


99.63% 


99.96% 


9. 


■nihil) & 1-0- i?ZZ;77 


0.07% 


0.00% (0 points) 


10. 


mh{2) & 1-0- Rzzrn 


0.00% (2 points) 


0.00% (0 points) 


11. 


/^77 & fJ-zz 


0.00% (2 points) 


0.01% (5 points) 



Table 2: Summary of the relative number of surviving posterior points, from the pre-LHC Bayesian global 
fits of the pMSSM, after imposing the Higgs discovery data. 

As we have already pointed out, the naive imposition of f-ijj (or, for that matter, any of the fix) 
as a constraint is fraught with danger in view of the uncertainties in the production cross section. On the 
other hand, the use of Rzz ;77 (see eq.©) largely eliminates these uncertainties, and given its sizable 
deviation from unity, is expected to lead to more robust claims. As Table|2]shows, the latter constraint is 
not significantly weaker than the former but, of course, it is far less constraining than the simultaneous 
imposition of both /x^^ and fizz- 

3.1 Implications 

The effects of the Higgs data on the pMSSM sparticle posterior distributions are shown in Fig. |2] which 
shows the allowed regions of sparticle masses for the model points that survive constraint no. 6 in Table|2] 
Note that the distributions for the surviving points are not meant to represent the possible outcome of a 
full-fledged parameter fit with the Higgs data, since this would require more statistic and convergence. 
Observations about the effect of the different constraints on the sparticle masses are as follows. 

• Only rrih = 122 — 128 GeV as a constraint: On imposing only this, about 40% of the posterior 
samples survive and the constraint does not affect the shape of the SUSY-breaking parameters and 
sparticles posterior distributions. The only exception is for the strong preference for multi-TeV 
trilinear coupling, l^^l ~ 5 TeV. The approximate conservation of the posterior distributions 
reflects the consistency of the pre-LHC fits with the Higgs boson mass. 

• Only Rzz u-y as a constraint: Although this discriminant shows an experimental central value 
very different from unity, note that it has yet a large error bar. Consequently, if we only demand 
a 2o" agreement with it, then essentially all the points survive for both log and flat prior samples. 
This is so because the pMSSM predictions are centred around unity, which is marginally consistent 
at 2o". The situation changes dramatically, if we strengthen the demand to a 1.67o" (con^esponding 
to a 90% C.L. for the 1-D marginalised distribution) agreement. Now, only about 13.3% of the 
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log prior and 0.4% of the flat prior posterior samples survive. Consequently, we skip the flat prior 
samples and continue the analysis with only the log prior. The surviving model points show a 
slight preference for lighter staus, stau-neutrino and lst/2nd generation squarks, and for slightly 
heavier ^4, > stops and 62- All other sparticle distributions remain approximately unchanged. 

• ruh = 122 - 128 GeV as well as Ezz^-y at 90%C.L.: 5.3% of the log-prior and 0.2% of the 
flat-prior posterior samples survive. The sparticles' posterior distributions for the initial log-prior 
sample set and for the surviving points are shown in FigiH As can be seen from the plots, there 
is a shift in probability towards heavier neutralinos (except for the LSP), charginos and stops. On 
the other hand, the data show preference for rather lighter lst/2nd generation squarks. The signal 
strengths' posterior distributions for the surviving points are shown in Fig. [3] 

• Only Rzz 177 at la: While a demand that a particular observable agree at la, statistically speak- 
ing, makes little sense, we, nonetheless examine this for curiosity. Partly, this exercise is motivated 
by the realization that, the error in, eq.©, has most probably been over-estimated, and is likely 
to go down significantly once the common errors have been factored out. Furthermore, with more 
data (already on tapes now), the statistical component would reduce substantially. However, work- 
ing with the currently available data, imposing this implies that only 0.3% of the log prior posterior 
samples survive while all of the flat prior samples get ruled out. The surviving points show prefer- 
ence for light (one to few hundred GeVl) staus and stau-neutrinos; preference for heavier t2> and 
moderate tan /3(< 40); and preference for heavy (> 500 GeV) X3 and (> 1000 GeV) X4 and xf- 
The mf_^ ^ posterior distribution remain approximately unchanged, mainly peaked around 2 TeV. 
With TUh = 122—128 GeV and Rzz-m 1"^' surviving points have the characteristics of: 
(/) having x? ^ "ih/2, having light (< 500 GeV) quasi-degenerate light lst/2nd-generation 
squarks with the LSP, and (///) light stau states. 

4. Summary and conclusions 

The signal strengths for the recently discovered Higgs boson-like state are used to propose a BSM dis- 
criminant, the ratio of diphoton- to ZZ- channel signal strengths. This has the advantage that sources of 
theoretical uncertainties cancel almost completely unlike the case for the individual signal strengths. The 
discriminant is applied to the posterior distributions for a pre-LHC Bayesian fit to the pMSSM, sorting 
for the impact of the discovery data on the sparticle mass predictions. The data provides mainly two 
classes of information; one about the mass and the other about the signal strength of the new state. 

The pre-LHC pMSSM fits predict an approximately prior-independent Higgs boson mass to be 
between 117 GeV and 129 GeV at 95% Bayesian confidence level and, hence, are in agreement with 
the recent discovery of a new state around 126 GeV. The pMSSM's signal strengths are naiTowly dis- 
tributed around unity, as shown in Fig. [T] suggesting compatibility with the SM expectations. These are 
in agreement with the experimentally determined values, to within 2a (or 95% C.L.) for the diphoton 
channel. This also holds for the aforementioned discriminant, simply as a consequence of the as yet 
large error bars associated with it. In other words, at this level, the effect of the signal strengths is not yet 
conclusive, and better precisions in the results is necessary for making distinctive conclusions. 

On the other hand, if we demand that the MSSM expectation of the discriminant agree with the 
data to within la, almost none, viz. 0.0% (3 model points) of the flat prior and some 0.3% (125 model 
points) of the log prior posterior samples survive. Indeed, even at 90% C.L., the constraints are quite 
strong. The surviving regions of the parameter space typically are associated with lighter stau states and 
lighter lst/2nd-generation squarks, whereas the stops, the charginos and the (non-LSP) neutralinos tend 
to be heavy. The prediction for light squarks is apparently in tension with limits from SUSY searches. 
However, it should be realized that the purpose here is to look for the directions in sparticle mass space to 
which the Higgs data pull, and the tension is not an overwhelming one. As for the stau states, the current 
direct LHC limits from either (Drell-Yan like) pair-production or staus from weak gaugino decays are 
weak. Future limits on the stau mass would be important constraints for the la Rzz ;77 surviving points. 
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Fig. 2: The plots compare the log prior pMSSM sparticle masses' marginalised 1-dimensional pre-LHC 
posterior distributions (dashed lines) and the surviving parameter regions (solid-blue lines) after impos- 
ing both ruh = 122.0 — 128.0 GeV and 1.67-cr Rzz i-y-y- All the masses are in TeV units. The vertical 
axes represent the relative probability weights of the model points. 
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Fig. 3: Marginalised posterior distributions for the signal strengths at 8 TeV LHC for the pMSSM points 
that survive the rrih = 122 — 128 GeV and 1.67-0" Rzz constraints. The colour spread from light- 
pink to blue represents relative probability spread from zero to unity. The inner and outer contours 
represent 68% and 95% probability regions respectively. The distributions for the flat prior sample is 
approximately the same but with narrower peaks around unity. 
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In conclusion, assuming an R-parity conserving MSSM, the results of our analyses imply that 
if the diphoton signal strength anomaly persists even after the accumulation of more LHC data, then 
it is very likely that there exist light sparticles which have to-date escaped detection. This could be 
due to the low magnitude of missing energy in the case of light (compared to the experimental missing 
transverse energy cuts) neutralinos and/or due to the quasi-degenerate nature of the spectrum which leads 
to soft jets Il24ll25l that would be buried under QCD background d Given the relatively low statistics 
for the surviving points, it is now time for updating statistically robust trends for MSSM physics by 
performing new pMSSM explorations and global fits to the post-discovery LHC data. Before this, short- 
range outlook indicates that performing a similar study, but for the impact of the LHC SUSY limits on 
the pMSSM posterior distributions is well within reach 1281 . The same is the case for checking the 
impact on the freedom of choice for considering the LSP as being either the only or a partial source for 
CDM relic density within the pMSSM framework |[29l . Other interesting directions include searching 
for plausible effects of non-diagonal flavour on the pMSSM predictions and for computing the impacts 
of the discovery data and SUSY limits on a phenomenological next-to-MSSM, the pNMSSM. 
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